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Introduction 


In order to obtain the data necessary to design future aeroassisted space transfer vehicles, 
NASA is currently developing the Aeroassist Flight Experiment (AFE). The vehicle used 
in this experiment will be deployed from the Shuttle Orbiter (scheduled for 1996), make a 
data-gathering aeropass through the upper atmosphere, and then return to orbit for pickup 
by the shuttle. The primary purpose of the AFE is to obtain data in the flight regime where 
chemical and thermal nonequilibrium effects dominate the shock layer, because such effects 
may play a large role in the design of aeroshells for future aeroassisted space transfer vehicles. 
A review of aeroassist concepts is given by Wa.lberg. 1 

The absence of applicable flight data and the inability of ground-based wind tunnels to sim- 
ulate the low-density, high-energy flight environment 2,3 is the impelling force which dictates 
the need for the AFE. This lack of data also creates an uncertainty in the heat transfer pre- 
dictions which are needed to design the AFE heatshield. Among the important factors which 
must be considered in the heating predictions are both chemical and thermal nonequilibrium, 
an extensive viscous region, small but finite wall recombination rates, nonequilibrium radia- 
tive heating, and three-dimensional effects due to a nonaxisymmetric heatshield geometry. 

Rochelle et al. 4,J have used an axisymmetric chemically reacting boundary layer code 6 to 
predict heating on the forebody heatshield of the AFE vehicle. However, a bounclarv-layer 
method does not properly account for the viscous-inviscid interaction and entropy- layer- 
swallowing effects. 

Hamilton et al. 7 used an axisymmetric chemically reacting viscous-shock- layer code 8 at lower 
altitudes and a Navier Stokes code 9 at higher altitudes to calculate the stagnation point 
heating rate for a typical ATE trajectory. This approach automatically accounts for viscous- 
inviscid interaction effects not accounted for in the above boundary layer results. The effect 
of wall slip was also investigated at high altitudes, but it was found to be unimportant for 
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the AFE. Although this approach, provides the best set of heating calculations for the AFE 
vehicle’s stagnation point it provides no information about heating for the remainder of the 
forebody heatshield. 

The purpose of the present paper is to calculate the total (convective and radiative) heating 
over the entire forebody heatshield of the AFE vehicle. The convective heating is calculated 
using a three-dimensional Navier-Stokes code (LAURA) which includes both chemical and 
thermal nonequilibrium effects. The output from the flowfield calculation is then used to 
provide inputs to a nonequilibrium air radiation code (NEQAIR) 10 to calculate the nonequi- 
librium radiative heating. Results are presented at two points on the current Baseline 5A 
trajectory corresponding to the start of the primary data taking period and peak heating 
(near the end of the primary data taking period). 


Symbols 


M 

Mach number 

P 

pressure, Njm 2 

q 

convective heating rate, W/cm 2 

s 

surface distance from nose along meridional plane, 

T 

static temperature, deg K 

Uoo 

free-stream velocity, km / sec 

x,y,z 

cartesian coordinates (see Fig. 1), m 

e 

emissivity 

0 

meridional angle (see Fig. 1), deg 

<7 

Stefan-Boltzmann constant, 5.6696xl0 -12 W/cm 2 K 

P 

density, kg/m 3 

Subscripts: 

c 

convective 

r 

radiative 

req 

radiative equilibrium 

total 

total 

w 

wall 

oo 

free stream 


Vehicle and Trajectory 


The aerobrake on the AFE vehicle (Fig. 1) is an elliptic cone, blunted with an ellipsoidal nose, 
and raked off at an angle of 17° relative to the cone axis. The ellipticity of the cone is such 
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as to produce a circular base in the rake plane which produces a highly three-dimensional 
stagnation region. A skirt is fitted tangent to the elliptic cone at the rake plane to reduce 
heating in this region. The aerobrake is made of a conventional aluminum structure covered 
with a thermal protection system with the outer layer consisting of tiles similar to those used 
on the Shuttle Orbiter. 

The vehicle is deployed from the shuttle payload bay. It is maneuvered to its entry atti- 
tude and then driven into the outer fringes of the atmosphere at a high velocity by a solid 
rocket motor which is jettisoned prior to the data-collection phase of the experiment. The 
vehicle enters the atmosphere at approximately 9.9 km/sec and penetrates to a perigee of 
approximately 75 km where its Velocity has decreased to 8.7 km/sec (see Fig. 2) It is then 
maneuvered back into a higher Earth orbit using aerodynamic lift and is retrieved by the 
Shuttle Orbiter and returned to Earth. The most important part of the mission, from a sci- 
entific view point, is the entry phase from 90 km down to 75 km where several instruments 
will take data. 

In the current paper, the heating is calculated at two points on the Baseline 5A trajectory: 
(1) peak heating and (2) the beginning of the quiescent period (see freestream flow conditions 
in Table I). The quiescent period is a period of approximately 30 sec when the attitude control 
jets are not fired (unless absolutely necessary to maintain control of the vehicle) so that data 
can be measured with minimum jet interference. 


Computational Methods 

In the present paper, convective heating to the AFE heatshield is computed using the 
LAURA code and radiative heating is computed using the NEQAIR code. A brief description 
of these codes is presented in this section. 

LAURA 

Gnoffo 11,12 has developed a three-dimensional, Navier-Stokes flowfield code, LAURA (Lang- 
ley Aerodynamic Upwind Relaxation Algorithm), which includes both thermal and chemical 
nonequilibrium. The LAURA code is described in detail in Refs. 11 and 12. The version of 
the LAURA code used in the present paper is a three-dimensional, second-order accurate, 
finite-volume algorithm that solves the thin-layer Navier-Stokes equations for laminar flow. 
The code uses a two-temperature thermal model and an 11-species chemical model for air. 11 
The thermodynamic and transport properties are taken from Ref. 14. 

The governing equations are relaxed in "pseudo" time until the solution reaches steady-state. 
The treatment of the governing equations is described as point-implicit because variables at 
the cell center of interest are treated implicitly, while the latest available data are used for 
the left-hand-side numerics. With this strategy, updating cell-centered variables requires the 
inversion of only a 16 x 16 matrix and only one level of storage. 
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The solutions presented in the present paper were computed with a grid of which had 64 
points between the body and the outer computational boundary. The grid on the body 
surface is shown in Fig. 3 where there were 42 points in the i direction and 21 points in the j 
direction. The surface grid was set up in this fashion instead of in meridional planes passing 
through the nose to avoid a “pole type” singularity at the nose. 

NEQAIR 


The NEQAIR 10 code is a nonequilibrium air radiation code that calculates the population 
distribution of atoms and molecules among their states by use of the quasi-steady-state 
assumption (QSS). The QSS assumes that the rates of population and depopulation of 
a given electronic state are much larger than the difference between the two rates, and 
the difference can be effectively ignored. Given the known population rates, the number 
densities of the excited and ground states can then be calculated from the predetermined 
thermodynamic state for any point in the flowfield. The NEQAIR code then calculates the 
intensity of radiative emission through a line-by-line integration over the spectral region in 
question. 

For the purposes of this investigation, the gas was assumed to be opaque below 2000 
angstroms and non-absorbing between 2000 and 15000 angstroms. Higher wavelengths were 
ignored. The total global emission was then integrated over the distance normal to the wall 
using the tangent-slab approximation to obtain the total radiative flux to the wall. 

The data that was input to NEQAIR (the species number density and temperature profiles) 
is obtained from the converged flowfield. 


Approach 


The forebody heatshield of the AFE vehicle is subjected to high heating on the entry phase 
of the mission prior to perigee. In order to design a vehicle that can survive and perform its 
mission successfully, one must accurately predict the heating environment. The conditions 
are such that the flow within the shock layer is in both thermal and chemical nonequilibrium. 
The surface heating is primarily convective, but the nonequilibrium radiative heating is too 
large to ignore completely. 

The LAURA code is used in the present paper to calculate the convective heating. Wall 
slip was neglected since it has been shown 7 that it is unimportant for these conditions. The 
surface of the aerobrake is covered with tile, made from the same material as that used on the 
shuttle orbiter, which has a small but finite chemical recombination rate. To account for this 
effect, the surface catalysis results of Stewart 15 have been used in the present calculations. 

Once a converged flowfield solution has been obtained, the species number density and 
temperature profiles at a grid point on the body are input to the NEQAIR code and the 
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radiative heating at that grid point is computed. This process is repeated at every second grid 
point on the body (in both directions) until the entire heatshield surface has been covered. 
The radiative heating at the intervening surface grid points is obtained by interpolation. 

The wall temperature used in the present analysis was the radiation equilibrium value (T„, t ,. e9 ) 
iteratively calculated from the following equation 

T w , req = [(<W + <w)/e<7] 1/4 (1) 


where the surface emissivity t was assumed to 0.85. 

The procedure for carrying out the computation was as follows: 


1. A initial wall temperature is assumed (usually constant or from a previous converged 
solution at different flow conditions). 

2. The LAURA code is run to obtain a flowfield solution and values for the convective 
heating rate (q^c). 

3. Using the species-number-density and temperature profiles extracted from the flowfield 
solution, NEQAIR is run to obtain the radiative heating (q™, r ) at each body grid point. 
This step need only be done once during the iteration process since surface temperature 
has a very weak affect on radiative heating using the present uncoupled approach. 

4. Equation (1) is now solved for a new value of wall temperature at each body grid point. 

5. A check is made for overall convergence by comparing the new wall temperatures with 
those used in step 2. 

6. If the wall temperature has not converged, go back to step 2 and repeat the process 
with the new wall temperature. 


The iteration process for wall temperature is stopped when the maximum difference between 
input (at step 2) and output (at step 4) wall temperature is approximately 50°K. This 
required three passes for the results presented in the present paper. 


Discussion of Results 


Distribution of surface heating along meridional planes for the peak heating point on the 
trajectory are presented in Fig. 4. The vehicle was at a = 0°. The flow conditions for this 
case is presented in Table I. On the lower portion of the body (<f> < 0° - Figs. 4(a)-4(f)) 
the heating decreases as away from the stagnation point (s ss 0), levels out over the conical 



portion of the body and then rises again as near the beginning of the skirt. The rise in the 
heating near the beginning of the skirt is caused by the rapid acceleration of the flow in this 
region. On the upper portion of the body ( <f> > 0°) the conical region decreases rapidly (see 
Fig. 1) and thus the tendency for the heating to level out on the conical segment of the body 
is less pronounced. The rise in heating at the beginning of the skirt is also less pronounced 
than it was on the lower portion of the body. 

The wall temperatures along meridional planes for this case are presented in Fig. 5. These 
are the resulting wall temperatures after three iterations as discussed in the previous section. 
The heating contours and wall temperature contours for this case are presented in Figs. 6 
and 7, respectively. 

Distribution of surface heating along meridional planes for the beginning of the quiescent 
pei iod are presented in Pig. 8. The vehicle was at a = 0°. The flow conditions for this 
case is presented in Table I. The heating follows similar trends to the peak heating results 
discussed previously but the overall level is somewhat lower because the freestream density 
is lower. 

The wall temperatures along meridional planes for this case are presented in Fig. 9. The 
heating contours and wall temperature contours for this case are presented in Figs. 10 and 
11, respectively. 

The surface heating rates have been interpolated onto the body point locations identified in 
the Aerothermodynamic Data Book 5 and these results are given in Tables II and III for the 
peak heating and beginning of the quiescent period, respectively. The results are provided 
in this form to assist with the heatshield design. 


Concluding Remarks 


The nonequilibrium convective and radiative heating rates on the heatshield of the AFE 
vehicle have been computed for two points (peak heating and the beginning of the quiescent 
period) on the baseline 5A trajectory. These results which are presented in both graphical 
and tabular form represent the most complete set of heating calculations available for the 
AFE vehicle. 
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Table I. Freestream flow conditions. 


Case 

a 

deg 

poo 

kg/m 3 

v 

y OO 

km/sec 

T 

OO 

°K 

Poo 

N/m 2 

Beginning of Quiescent Period 

0 

2.176xl0 -5 

9.715 

196.8 

1.227 

Peak Heating 

0 

3.741xl0~ 5 

9.306 

198.9 

2.134 
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Table II. - Heating rates on body at peak heating. 
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Table II. - Concluded. 
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Table III. - Heating rates on body at beginning of quiescent period. 
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Figure 1 . - AFE vehicle geometry. 
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V, km/sec 

Figure 2. - Baseline 5 A trajectory. 
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Figure 3. - Body Grid. 
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(b) <J) = -80° 
Figure 4. - Continued. 
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Figure 4. - Continued. 
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(d) <|> = -40° 
Figure 4. - Continued. 




3 W o/M ,M b 


22 


(e) <|> = -20° 
Figure 4. - Continued. 
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(f) <t> = 0° 

Figure 4. - Continued. 
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Figure 4. - Continued. 
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(h) 0 = +40° 
Figure 4. - Continued. 




2 ujo/m ,M b 


26 


Figure 4. - Continued. 
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G) <t> = +80° 
Figure 4. - Continued. 
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(k) <\> = +90° 
Figure 4. - Concluded. 
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Figure 5. - Continued. 
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(c) <j) = -60° 

Figure 5. - Continued. 
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Figure 5. - Continued. 
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(e) <j> = -20° 
Figure 5. - Continued. 
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Figure 5. 
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(k) <j) = +90° 
Figure 5. - Concluded. 
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Figure 6. - Heat Transfer Contours at Peak Heating Point. 
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Heating Point. 



2 iuo/m ,M b 


44 


Figure 8. - Heat Transfer at Start of Quiesce 
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(b) (j> = -80° 
Figure 8. - Continued. 
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(c) <j> = -60° 

Figure 8. - Continued. 
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(d) <\> = -40° 
Figure 8. - Continued. 
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(e) <t> = -20° 
Figure 8. - Continued. 
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(f) 4> = 0° 

Figure 8. - Continued. 
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Figure 8. - Continued. 
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Figure 8. - Concluded. 
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Figure 9. - Continued. 
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(c) <j) = -60° 
Figure 9. - Continued. 
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Figure 9. - Continued. 
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(e) <)> = -20° 
Figure 9. - Continued. 
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Figure 9. - Continued. 
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Figure 9. - Continued. 
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Figure 9. - Concluded. 
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of Quiescent Period. 
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